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I. INTRODUCTION 

Hadron spectroscopy in lattice QCD is entering a new 
era, in particular, recent developments in the application 
of variational methods [IHS] to large bases of hadron in- 
terpolating fields have made the extraction of the excited 
spectrum of hadronic states a realistic possibility (see e.g. 
[IHS]). Since excited hadrons appear as resonances in the 
continuous distribution of multi-hadron scattering states, 
to study hadron spectroscopy one requires evaluation of 
scattering amplitudes, but because lattice QCD is for- 
mulated in Euclidean space, we do not have direct access 
to these [7] . Fortunately, in a finite volume, interactions 
between particles as they evolve from the in to the out 
states lead to discrete changes in a free particle's energy 
that can be related to the scattering amplitude [8j. 

Various extensions to the framework derived by 
Liischer in ^ have been proposed which allow for evalu- 
ation outside the center-of-mass frame [51-fT5]. and to in- 
clude the coupled-channel effects that can appear above 
the inelastic threshold [Tmi8j . The original approach 
and its extensions to describe the moving center-of-mass 
frame have been quite successfully used by the lattice 
community to extract elastic hadron-hadron scattering 
phase shifts [S [IMS]- 

In this work, we discuss a generalization of Liischer's 
method for relativistic scattering in terms of a Hamilto- 
nian where the specific interactions considered are based 
on a relativistic particle exchange model. We apply 
the technique to a two-channel system and a general- 
ized Liischer's equation for scattering in a moving frame 
is derived based on the relativistic Lippmann-Schwinger 
equation. The coupled- channel system has been consid- 
ered previously, [TH HSl [IZl US], and our result agrees 



with the form derived in a non-relativistic formalism in 
[Tl] . A novelty of the present work is to discuss practical 
strategies for extraction of scattering amplitude param- 
eters (phase shifts and inelasticities) from lattice sim- 
ulations of a coupled-channel system. These strategies 
are demonstrated using an explicit toy model of resonant 
two-channel scattering. 

We also discuss the three-particle system, consider- 
ing the finite-volume representation of the isobar model 
[?7] in which interactions in the three-particle sys- 
tem are approximated by two-body scattering. We de- 
rive expressions which relate finite-volume energy shifts 
to isobar-model scattering amplitudes. 

The paper is organized as follows. A discussion of 
the elastic scattering in a finite-volume is given in Sec- 
tion Inl with the extension to the coupled channel sys- 



tem in Section HI Strategies for extracting scattering 
amplitudes from measured discrete finite- volume spectra 
are presented in Section |TVj In Section |V] we discuss the 
three-particle system, with summary and outlook given 
in Section IVll 



II. 



FINITE- VOLUME ELASTIC SCATTERING 
IN A HAMILTONIAN FRAMEWORK 



In this section we present relativistic two-particle scat- 
tering on a torus using the Hamiltonian formalism de- 
veloped in |281 129j . In particular we consider a complex 
scalar field, describing a charged boson, (p^, of mass 
m, and its interactions with a neutral boson, 0, which 
acts as a force carrier and is described by a real scalar 
field, 8. We first derive the Liischer formula describ- 
ing the finite-volume spectrum of the asymptotic two- 
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particle, 



det 



state, 



5jM,j'M' cot5j{k) ~ Mf^jj,M,{k) 



0, 



where the volume and scattering-frame dependent ma- 
trix element A^j*^) ji^, is defined in Eq.(Bl) and (B3|, 
and the center-of-mass frame scattering moment urn, k, 



is related to the energy by - 
responds to a Lagrangian density. 



i^. The model cor- 



. (1) 

from which the Hamiltonian can be derived following the 
canonical procedure of instant-time quantization (see Ap- 
pendix |A]) [30] . Taking matrix elements of the Hamilto- 
nian in an infinite basis of Fock states spanned by any 
number of 4> and 9 bosons one can obtain a Schrddinger 
equation H\'^) — E\'i/) for the eigenstates of the theory. 
Assuming /x 3> m, in describing low-energy (/)-boson scat- 
tering we can truncate the Fock space to include up to 
one 0-boson in the intermediate state, which reduces the 
Schrodinger equation to 



H22 H23 
H32 H33 



= E 



\cjy+rO) 



(2) 



The three-particle sector can be formally eliminated, re- 
sulting in an effective two-body equation, 



[E - i?22)i0+r ) = H2. ^„ H32W 



(3) 



A. Two-particle scattering in an infinite volume 

Before considering two-particle states on a torus, we 
will first review the scattering problem in infinite- volume, 
with further details given in Appendix [X] After elimi- 
nating the three-particle states |(/)+(?5>^0) from the cou- 
pled system (c/Eq.([2])) we are left with an equation for 
the center-of-mass frame momentum-space wavefunction, 
(/jjM(q), which is a product of a radial wavefunction de- 
pending on the magnitude of the relative 3-momentum, 
q = |q|, and the spherical harmonic of definite angular 
momentum, {J,M), 



'/fjA/(q) 



1 



2v^ 



(4) 

Here, E = y/s is the energy of the two-particle system 
in the center-of-mass frame. The non-local potential, 
T^(q, k), induced by 0-exchange is given explicitly in Eq. 
(Al). Expressing this equation in coordinate space via a 



Fourier transform gives 



JM 



(r)- j(fv'Go{v-v'-^s) jd\V{v',-^)i,jM 



(5) 



where the free Green's function is given by 

^3q giq(r-r') 



Go(r-r';Vi) = 



(27r)3 y^-2yq5 



(6) 



The wavefunction satisfies a relativistic Schrodinger 
equation. 



(7) 

While Eq.Q was derived in the context of a particular 
model, our subsequent derivation only requires the gen- 
eral form of the relativistic Lippmann-Schwinger equa- 
tion. The asymptotic component of the two-body wave- 
function relevant to scattering is given by the large dis- 
tance behavior of the Green's function. Evaluating the 
integral in Eq.([6]) (c/. Appendix [A| , we find 



Go(r; 



v/i-2v/A^ 
e^'^'' 1 
2 Anr r 



' Pdp 
(2^)2 



-pr 



fc2 + p2 



(8) 



with the first term on the right hand side dominating 
as r -> cx). For a potential V which falls at large sepa- 
rations, the solution to Eq Q outside the range of the 
interaction is given by 



V'JM(r) ^ l^i'YjMi^) [^^j.,{kr)+ikfj{k)h+{kr)\ 



(9) 



where f.j{k) is the partial wave scattering amplitude. 



- -jjJdh'd'zj,j{kr')Yj\j{r')V{r',~z) Vja,(z), 

(10) 

which up to the inelastic threshold can be parameterized 
in terms of a single real momentum-dependent parame- 
ter, the scattering phase-shift, Sj{k), as 



fjik) 



4tt 



e'^-' sin 5, 



B. Two-particle scattering on a torus 

Now we consider the theory in a cubic box of volume 
V — , with periodic boundary conditions. In Eq. ^ 
we split the integral over r' into a sum of integrals over 
a set of boxes labelled by the integers n representing the 
location of one of its corners, 

X jd^z' V{r' + nL, -z' - nL) i}%{z' + nL). 

(11) 
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In general we can make the wavefunctions periodic up to 
a phase, 



nL) 



'/'jm(z)> 



where the Bloch wave-vector, Q, is related to the total 
momentum of the two- particle system [9] by P = 27Q. 
7 = Vs + P^/-\/s is the Lorentz contraction factor that 



reduces the effective size of the box in the direction 
parallel to P. Using the periodicity of the potential, 
V{y' + nL, — z' — nL) — V^(r',— z'), and the boundary 
condition on the wavefunction, we have 



(r)= / dVGQ(r-r';^/i) 

JL3 



which is analogous to the infinite-volume equation, but 
with the Green's function given by 

GQ(r - r'; ^/s) = ^ Go(r - r' - nL; ^s) e^Q "^ 

n6Z3 

Using the Poisson summation formula, 
obtain 



GQ(r-r';yi) = - ^ 



„iq-(r-r') 



y/s I ^ e«q-(r-r') 



where Pq = {q g M^] q = ^n + Q, for n e Z^}, and 
where we have retained only the leading term in the limit 
L^|r — r'l^rn"^. Finally, expanding Eq. (Ill for r ^ 
r' and using the definition of the scattering amplitude, 
Eq. ( 10 1 we can express the wavefunction as 



2m 



J'M' 



5jM,.j'M' nj,{kr) - M^jfj j,j,j,{k) jj-{kr) 



(12) 



The residual sum over all angular momenta reflects the 
broken rotational invariance induced by the finite cu- 
bic volume, with the volume-dependent matrix elements 
M. given in Appendix [B] In the infinite-volume case, 
the most general solution of the relativistic Schrodinger 
equation, Eq.([7]), outside the range of the potential is 
LjA/CjM V'jMi") for V'JM(r) given by Eq. (l9|. Cor- 
respondingly in finite- volume, the most generalsolution 

is given by J^jM'^JM^fi^Hr) for V'jM*^'(r) given by 
Eq.(12|. Matching the two wavefunctions at a fixed r. 



larger than the range of the interaction, we obtain 
Y,cjMYjM{r)f^ [4^jj(fcr) +ifc/j(fcr) h+{kr)] 



JM 



J2 C.,Ml'kJj{k)Yj,M'{v) 



JM,J'M' 

X \SjM,J'M' n,j>{kr) 



M 



(Q) 

JM.J'M 



(fc) jj, (kr) 



which has a non-trivial, cjm 7^ 0, solution provided 



det 



SjMj'M' cotSj{k) - M 



'-^^ (k) 



0. (13) 



This condition expresses the relationship between the 
asymptotic behavior of the two-particle wavefunction 
on a torus, expressed through the matrix elements 
M 

M^i,j'M' ' ^^'^ effect of the interaction on the wave- 
function determined by the phase shifts, Sj. In practice 
for a given set of elastic phase-shifts, 6j{k), it determines 
a discrete spectrum of states in a finite volume. 

The analysis presented here can be generalized to an 
arbitrary shaped box. In general the three edges of box 
are spanned by three arbitrary vectors Li 2.3. The vol- 
ume of the cube is replaced by (Li x L2) • L3 and the 
vector nL by 



2,3 ' 



UiTiijUi £ 7L. Finally the momen- 
tum q = 27rn/L, n e is replaced by generalized mo- 
mentum 27rX]i=i,2,3'^i(Lj X Lfc)/|(Li x L2) ■L3|,'«z G 
where indices (z, j, fc) follow the cyclic permutation. Such 
a generalization has to be considered when using the 
moving center-of-mass frame since the symmetric shape 
of a cubic box in the rest frame is deformed due to 
Lorentz contraction [9 . In this case if P = 27rd/L, d e 
Z^, is the center-of-mass momentum, the volume of the 
box becomes 7L'^, and the vectors nL and 27rn/L are 
replaced by 7nL and 27r7~^n/L, respectively (using the 
notation defined in P|). With these substitutions and 
the relation P = 27Q, our definition of the matrix el- 
ements A^7*^| //jv// becomes identical to the matrix ele- 
ments Mf^^,'^, in Eq.(89) of 1. 

Typically, as discussed in [8j , for the low-energy region 
that we are interested in, higher partial waves become 
progressively smaller and can be ignored, so that the par- 
tial wave basis can be truncated at a certain maximal 
angular momentum Jmax- For a finite- volume with cubic 
boundaries, the continuous rotation symmetry is reduced 
to the little group of allowed cubic rotations that leave 
the centre-of-mass momentum invariant - the matrices 



in Eq.(13l become block-diagonal if subduced according 



to the irreducible representations of these little groups. 
Details of subduction in general moving frames can be 
found in [311. 



III. COUPLED CHANNEL SCATTERING IN 
FINITE VOLUME 

We extend the model of the previous section to include 
additional two-particle asymptotic states, by adding an- 
other species of charged bosons, cr*, which also couples 



4 



to the force carrier, 6*, into the Lagrangian. We can ob- 
tain coupled equations for the two-particle states, |(/>^0^) 
and Ict^ct^) , by eliminating states featuring three or more 
particles and obtain a two-channel Schrodinger equation, 



1^+0-) 
\a+<j-) 



1 



E-H. 
1 

E-H, 



[y.^|0+0-)-f y..|c7+a-)] , (14) 



where H^^\ H^"^ are the one-particle operators and 
V^^,V^cr,V„^,Vcra are effective interactions (two-body 
operators) generated by the reduction to the two-particle 
subspace. From Eq.(14|, for the channel wave functions, 



V'jM = (''1'^' •^-^)' obtain 

rjMir)^ y'dVG^(r-r';v^) 

The coupled-channel scattering amplitudes can be de- 
fined by 



m 
i 



f Jdh'd^zjj{k^r')YJ,i{r') 

xV^p{r'-z)^^j^,{z), (15) 



where ka = •\/s/4 — is the magnitude of the relative 
momentum in the center-of-mass frame of the two parti- 
cles in channel a. By analogy to the single channel case. 



the asymptotic wavefunction in channel a is given by 



2m, 



-YjM{r)i-' 



^7rjj(k^r)+tk^hj{k^r)j:pff{s) . (16) 



The two-channel scattering matrix f"^ is conventionally 
parameterizated in terms of two scattering phase-shifts, 
Sj{s) {a = (/>, (t), and an inelasticity, r]j{s), representing 
the fraction of flux exchanged between the two channels. 



/r(5) 



Att rjj e 

kfy 



2iS°, 



2i 



ff{s) = 



Extending the one-channel analysis of the asymptotic 
states in finite-volume to the two-channel system, one 
obtains. 



JM 



J'M' 

X SjMj'M' nj'ikar) - M'jfjj,]^j,{ka) jj'ik^i 

(17) 

Matching the wavefunctions in finite- volume, Eq. (fT7|, 



to the wavefunctions in infinite volume, Eq. ( |16[ ), we can 
derive the analogue of Eq. ( 13 ) for the coupled channel 
case. 



det 



JJM.J'M' 



cot A* - M 



(Q) 

JM.J'M 



iS. 



JM,J'M' 



M 



(Q) 

JM,J'M 



2r(,7 sin AJ 



« 0JM,.J'M' 



M 



(Q) 



1M,.TM' 



KIM, .I'M' 



cot A} 



M 



(Q) 

JM,J'M 



2r]j sin A*^ 



= 0, (18) 



where A}(s) = S'}{s) 



One can show that 



^ log rj J (s). 

this result is equivalent to Eq.(4.14) in [14 . As in the 
single-channel case, the subduction of Eq.(18l to irre- 



ducible representations of the appropriate little groups 
can be performed [5T] . 

Eq.(18) is complex and leads to two equations for the 



three real parameters, Sj,Sj,rij, in each partial wave. 
Further, under certain conditions, e.g. S-wave domi- 



nance, Eq.( 18 1 becomes a single real equation in the three 
unknowns, ((5q , (5g , 770). Thus additional constraints need 
to be imposed to obtain a unique solution. For example, 
in [15] unitarized chiral perturbation was used to con- 
strain amplitude parameters at low energies. In the next 
section we explore other strategies in the context of an 
analytical parametrization of the amplitude. 



IV. A TOY MODEL OF RESONANT 
COUPLED-CHANNEL SCATTERING IN 
FINITE-VOLUME 



In order to explore possible strategies for extracting 
coupled-channel scattering amplitudes from the discrete 
finite-volume spectra emerging from lattice QCD com- 
putations, we consider a simple model of two-channel 
S-wave scattering. The model is based on resonance- 
dominated scattering and satisfies the analytical proper- 
ties required of partial wave amplitudes. With an ex- 
plicit model for two phase shifts and the inelasticity we 
can solve Eq. ( 18 ) to obtain finite- volume spectra of 
states as a function of the volume {V = L^) and total 
momentum of the center-of-mass P = 27rd/L, d € 
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FIG. 1: Phase-shifts and inelasticity for the model defined in 
the text. 



(P = 27Q). We then use this spectrum as pseudo-data 
representing a hypothetical lattice QCD simulation and 
attempt to reproduce the input model. 



A. Analytical model of two-channel scattering 

We consider a model in which a single S'-wave reso- 
nance in both channels interferes with a non-resonant 
background. The two-channel scattering amplitude is 
parametrized in terms of a if-matrix, 



gag 13 



(0) 



M2 - S 

which is related to the t-matrix by 



(19) 



and to the scattering amplitude defined in Eq. ( 15 ) by 

ij'^i^) = V^fj^{^)/^^ ^ith = 0, o"). Here Ia{s) is 
the Chew-Mandelstam form [32], whose imaginary part 
above threshold, i.e. for s > 4m^, is given by the phase- 
space, 



= /a(0) - - / ds'xil 



4m2 



1 



s' (s' — s)s' 



We have opted to subtract the integral once, and it is 
convenient to choose /q(0) such that Ke laiM"^) = so 
that we have an amplitude which for real s near is 
close to the Breit-Wigner form with mass M. The t- 
matrix is an analytical function in the complex s-plane 
with the discontinuity across the right-hand cut deter- 
mined by unitarity. 

With the following choice of parameters, 

= 0.2 GeV, = 0.4 GeV 
M = 1.1 GeV, 30 = 0.35 GeV, = 0.2 GeV 

^ 7?:^ = 0, 7i°j = 0.7, 7i^j = 0.7GeV-^ 7^"/^' - 0, 



we obtain the phase-shifts and inelasticity shown in 
Fig. [1] The parameters have been chosen in such a way 
that there is a narrow resonance near s = (1.1 GeV) 2. It 
is usual to analyse scattering in terms of the most rele- 
vant singularities of the i-matrix on the nearby unphys- 
ical sheets. Poles on unphysical sheets are often identi- 
fied with hadron resonances. In this case the four sheets 
(sheet I is the physical sheet) can be defined by 



[*sheet( 



2i 
2i 
2i 



4mj 



Sad 



4mi 



4nC 



sheet I 
sheet II 
sheet I II 



Saa sheet IV 



The model amplitude has a single pole on each of sheets 
II and ill, with the i-matrix in the neighbourhood of the 
pole at So behaving like, 



So - S 



with 



^ = (1.1072 - |0.0888) GeV 

C0 = (0.3511 GeV) e-* ™!^'^, c„ = (0.1339 GeV) g-'^ ^^S'^ 
on sheet II and 

0.1096) GeV 

= (0.1396 GeV) e+'0-249- 



(1.1083 - I 
= (0.3499 GeV) e+'" °°°9^ 



on sheet III. Our aim is to use the finite- volume spectrum 
determined on a set of volumes and total momenta, P, to 
reproduce the pole positions of this scattering amplitude. 



B. Finite- volume spectrum 

The finite-volume spectrum corresponding to the 
model defined in th e previous section can be obtained 
by solving Eq. ( 13 1 for energies E = y's < 2m^ and 



Eq. (18) for energies above 2mo-- Restriction to S'-wave 



scattering reduces Eq. ( [T8| to, 

= V[M^- M^] sin (^* - S^) 

-[l-M^M^lcos + 

where Mct> = -^/^o a/=o /'=o j\/'=o(^0) '^ith ^ similar ex- 
pression for Ma- Q is a function of d, L, E as dis- 
cussed in Sec. [TTbI In Fig[2] we show the finite-volume 
spectrum obtained by solving Eq. ( 13 ) and ( 18 1 as a 
function of the volume in a region L = 16 — 24 GeV^^ 
(or L = 3.2-4.7 fm). 



(20) 



FIG. 2: Finite-volume spectra for the A"-matrix model described in Sec. |IV A| Black dots indicate the spectrum obtained 
by solving Eqs. (13 1, (181. Red and blue curves represent the energy of a non-interacting pair of mesons (a — 

2 



(^^ml+kl + y/ml + ki) " 



1/2 



ki + k2 = P and k = f-n, n e 



C. "Pointwise" estimation of scattering from 
finite-volume spectrum 

One approach to determining the phase-shifts and in- 
elasticities at discrete values of energy is to locate mul- 
tiple energy levels (in different volumes and/or different 
d) that appear at approximately the same energy. As 
an example consider the three levels labeled A in Fig. [2] 
which all lie within 2 MeV of = 830 MeV. For the 
three levels we can build three independent copies of Eq. 



( 18 1 which each feature approximately the same values of 



d^s), (5°"(s), 77(3), which can be determined by solving the 
set of simultaneous equations. Since the energies are not 
exactly degenerate, there need not be an exact solution 
to the equations and hence we seek to find the solution 
which minimizes 



E 

E(L,d) 



levels near 1.117 GeV labeled d in Fig. [2j Linear inter- 
polation between the energies of the E and F points in 
Fig. § gives (5-^(1.117 GeV) = 102.3°. Using this value 
and minimizing with respect to (5°^, 77 at the energy cor- 
responding to point d we obtain 6°^ = 46.2°, 77 — 0.800. 
A spline interpolation using all the points A — G yields 
(5'^(1.117GeV) = 102.2° which results in = 46.3°, 77 = 
0.799 at the point d. In Fig. [3] we show the results for 
sets of two degenerate levels labeled a — f from Fig. [2] 
In each case the errorbars represent limiting values ob- 
tained using two methods of interpolation. Even though 
in some cases there is a considerable sensitivity to the 
interpolation method, overall the points are in reason- 
able agreement with the model input (solid light-colored 
curves) . 



D. Parameterized estimation of scattering from 
finite-volume spectrum 



with defined in Eq. 
the three energy levels. 



20|) and where the sum is over 
For the levels A, the obtained 



considered, E 
three volumes, 



0.1 

L 



solution, as shown in Fig. [sj is within 3% of the exact 
value of 6'*'{sa),S'^{sa),v{sa)- 
Within the energy region 
1.2 GeV, considering only 
16, 20, 24GeV~^, and three sets of center-of-mass mo- 
mentum, d = (000), (100), (110), we can isolate a num- 
ber of such sets. These sets of three near-degenerate 
energy levels are labelled A — G in Fig. [2j In Fig. [3] the 
labels are shown on the plot of 6'^ with the corresponding 
solutions for 5" and 77 marked by solid dots. 

With these points alone, in Fig. |3]we see strong hints 
of a signal for resonant behavior in the 6'^ phase shift. 
While obviously resonably successful, this approach does 
not make optimal use of the finite-volume spectral in- 
formation, by failing to use any energy level which does 
not have two near-degenerate partners. To use more of 
the discrete levels we might consider building a system of 
two instances of Eq. ( 18 1 with one parameter from the 
set 6'^,S'^,i], estimated using interpolation between al- 
ready determined values. For example consider the two 



The previously discussed "pointwise" strategy, 
while having the advantage of being largely model- 
independent, is reliant upon there being multiple energy 
levels which, through accident or design, are close to 
degenerate. Since it would be unusual to engineer lattice 
volumes purely for this purpose, it is appropriate to 
consider alternative methods of analysis. One such 
approach that makes full use of all determined levels 
involves parameterising the scattering amplitude and 
performing a minimisation to get the best description 
of the determined finite-volume spectrum by varying 
the parameters. In the current toy-model, even limited 
"pointwise" analysis would suggest that the phase S'^ is 
rapidly rising and would indicate that a resonance could 
be present. By including a pole (as well as polynomial 
behaviour) in a ii'-matrix we are likely to get rapid 
convergence to a solution with a pole in the i-matrix. 

We will take this opportunity to make the toy model 
a slightly more realistic simulation of an actual lattice 
QCD calculation by introducing statistical uncertainty 
on the energy level values. In recent work [52], [5], 
the Hadron Spectrum Collaboration has obtained sta- 
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FIG. 3: "Pointwise" determination of the phase-shifts and inelasticity. Energies A — G determined from constrained three-level 
analysis, energies a — f by interpolation in 5*" from two-level analysis. The light-colored curves show the exact input model. 



tistical errors on excited levels as small as 0.3% and we 
will assume that this remains practical. For each energy 
level below 1.2 GeV on volumes L = 16,20GeV"^ with 
d = (000), (100), (110), we randomly generated an en- 
semble by drawing from a distribution whose mean is 
the exact value given in Fig. [2] and whose variance is 
chosen such that the ensemble has variance on the mean 
of 0.3% of the mean value. The resulting spectrum is 
shown in Fig. |4] 



Parametrizing according to the form given in Eq. ( 19 1, 
we can minimise a function. 



X^({«J) = 



E 

B„(i,d) 



(where E'^°^ are the solutions of Eqs. (13) and (18l), 
by varying the parameters of the if- matrix, {ai\ — 
{M, f;(j,7(") . . .}. In practical lattice QCD calcula- 
tions, the can be trivially adjusted to deal with cor- 
related data by replacing the diagonal variance by the 
inverse of the data covariance matrix. 

In FigjS] we show the parameterized phase-shifts and 
inelasticity obtained using three models: 



A: "exact model", which uses a 1st order polyno- 
mial in s to describe the non-pole contribution to 
the i^T-matrix with 7^0"'^'' — 70^'^'' = [5 parame- 
ters] 



• B: "relaxed model", with 1st order polynomial in 
all channels i.e. all 7(0'^' free [9 parameters] 

• C: "tight model", with 0th order of polynomials in 
all channels i.e. 7^^-' = [6 parameters]. 

As one would expect, within statistical uncertainty, pa- 
rameterization A reproduces the input model. Parame- 
terisation B, which is more flexible, also reproduces the 
input quite well over the energy region where data is 
given, but begins to deviate from the original JiT-matrix 
in the energy range outside of the fit region. Parameter- 
ization C does not have sufhcient flexibility to describe 
the complete energy dependence, and it fails to describe 
S"^ at high energies, away from the energy region where 
the pole dominates. 

Our principal interest lies in identifying resonances as 
poles in the complex-s plane - analytically continuing the 
fitted model amplitudes we find that all three have single 
poles on sheets II and III whose locations are in a rather 
good agreement with the input pole position (see Fig. |6|. 
The residues at the pole agree similarly. 

In summary, the "pointwise" strategy may provide 
a less mo del- dependent approach for extracting phase 
shifts and inelasticities, however, this method is limited 
by the number of points for which accidental degenera- 
cies appear. Parameterizing the scattering amplitude al- 
low us to make use of all measured energy levels, how- 
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> 



cause of the low mass of pions, three-particle states be- 
come important at relatively low energies. Methods for 
extracting three-particle amplitudes from lattice simula- 
tions are thus also needed. In general, the scattering of a 
multiple-particle system is quite complicated, however, in 
some cases the quasi-two body approximation has proven 
to be quite successful 1331 131] ■ A recent consideration of 
the three-particle system on a torus in the framework of 
Faddeev equations is presented in |35j . 

In the presence of two-body resonances, the three- 
particle system can be approximately described in terms 
of isobars [261 EH EH [37] , for example, the decay process 
of i^i(1400) to the KiTTT final state can be quite well ap- 
proximated by /ir*(892)7r in S'-wave, with the K* isobar 
decaying to Kn in a P-wave. 

In this section, we will generalize the previously intro- 
duced methods to the three-particle system in the isobar 
approximation. For simplicity we consider three scalars 
with equal masses, m. Under the assumption of strong 
pair-wise interactions between two particles, the Hamil- 
tonian of the three-particle system takes the form 



H = 



fc=l,2,3 



pi + m? 



fc=l,2,3 



L= 16GeV" 



L = 20GeV" 



FIG. 4: Orange rectangles: finite-volume spectra with 0.3% 
noise. Black lines: exact finite-volume spectrum given in 
Figi 



ever we need to find suitable parameterizations. One 
strategy is to explore the "pointwise" approach to find a 
crude guide to the energy dependence and then build pa- 
rameterizations which are able to reproduce the obtained 
form. The parameterizations, which should respect cer- 
tain constraints applicable to scattering amplitudes, can 
be made progressively more sophisticated in an effort to 
reduce the overall - in this sense the approach is not 
dissimilar to what is done with real experimental data. 



V. THREE-PARTICLE SCATTERING ON A 
TORUS WITHIN AN ISOBAR APPROXIMATION 



where (p^, x^) are the momentum and coordinate of the 
fc-th particle. In the center-of-mass frame, the wavefunc- 
tion of the three-particle system has the following form 



V'Jm(xi,X2,X3) = ^ (j)[jj^j{rij,rk), 

fe=l,2,3 



where r^- = — Xj is the relative position between two 
particles forming the (zj) isobar and = ^(xi+Xj)— 
is the relative position of the spectator particle and 
the isobar. The solution of the Schrodinger equation 
HipjM = EipjM is given by 



X V*j(rjj)V'JM(x'i,x'2,xj5), 



In the previous section we considered only two-particle 
Fock states in the inelastic region, however in QCD, be- 



with Go being the relativistic three-body Green's func- 
tion. 
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FIG. 6; Poles of the f-matrix on sheets 11,111 determined from 
parameterizations A — C as described in the text. Black dots 
show the exact poles of the input model. 



(27r)3 (27r)3 



pjqij-rij+ikfc-rfc 



E- yfn?' 



+ (ikfc + q, 



In the isobar approximation, interactions within the 
three-particle system can be described by two types of 
process: interactions between the isobar pair and the 
spectator, [if) + fc o iij) + fc, with all possible arrange- 
ments, and interactions between particles within the iso- 
bar, i ^ j. This enables us to represent the full three- 
body Green's function by a product of two independent 
two-body Green's functions 



Go(i"jj,rfe;£') 



2n 



{vu-,E) Gl+\r,,;E,,), 



where 



G, 



{ij)+k 





{rk;E) 



d^k. 



(2^)' E - ^m^ + ^ - ^Ef^ + 



is the free Green's function for propagation of the system, 
of total energy E, made of the isobar pair (ij) (with 



invariant mass Eij) plus the spectator particle, and where 



d^qij 
(27r)3 



Ei. 



is the free Green's function describing the propagation of 
particles, inside the isobar. 

The Green's functions, excluding any pieces which de- 
cay exponentially with distance, are explicitly given by 



Go{rij,rk;E) 



1 
27ri 



E- 



dEij 



E (47r)2 
e 



(21) 
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with 



^[E^ - (m - £;,,)2] [iJ2 _ + i5^^.)2] 



2£; 



In Eq.(21), J^^dEij has been replaced by dEij, 



becauseThe three-body Green's function has oscillatory 
behavior for Eij only inside the physical region. The 
asymptotic behavior of the three-particle wavefunction 
is given by 



■!/'Jm(xi,X2,X3) 

- T 



fe = l,2,3 



dE, 

2m 



^ fjil-.S^jMs- ■ iLkMi,. ilE„ , kE.j), 

where the three-body scattering amplitudes are defined 

by 



xi-^'=jL,(fc£.,r^)r£^M,,(ffc) 

X V^,,(r^)^jM(x'i,x^,x;j). (22) 



The total angular momentum of the three-particle sys- 
tem is J, made up of (S'ij, Mg-^.), the spin and its z-axis 
projection of the isobar pair (ij), and {Lk,Mi^^), the rel- 
ative orbital angular momentum and projection between 
fc-th particle and isobar pair (ij). 

The partial wave amplitude can be parametrized by 



= (5., Ms,, ; LfeMi J JM) /s., {qe., ) /5.,l. j ), 

(23) 



where 



An 



fs,,{qE,,) = ^^e"^^'J- sin (5s 



is the partial wave scattering amplitude of two spinless 
particles inside the isobar pair (ij) and 



fs,,L,j{kE,,) = -^e^''^^'''''' sm6s,^L,j 

KEh 



is the partial wave scattering amplitude of the isobar pair 
(ij) and the spectator particle. 



The final expression for the three-particle wave function, including the homogeneous term, reads 



V'Jm(Xi,X2,X3) 



fc=l,2,3 SijMs. 



{S,,Ms^r,LkMLjJM)Ys^ 



E-ra 
2m 



dEij i- 



(47r) js,, (fe,, r,j) + i qE,^ h+ (qE,^ nj) fs,^ {qE,^ ) 



X i^" [(4^)jL,(fcB,,.rfe) h+^{kE^,rk) fs,,L,j{kE^,)\ 



When considering three particles in a cubic finite-volume we may choose the following boundary condition on the 
wavefunction 



where k = 1,2,3 and Q is the Bloch wave vector for the three-particle system. Boosting the three-particle system 
from the center-of-mass frame to a lab frame with the center-of-mass of isobar pair {ij) fixed at the origin, Xj -|- Xj — 
0, the wavefunction of the three-particle system at the lab frame can be written as the product of a plane-wave, 
giP (xi+x2+x3)/3^ g^j-^jj piece depending only on relative coordinates, Yij = x^ — Xj and = — x^. As in the case of 
two-particle scattering [9], requiring periodicity of the lab frame wavefunction with respect to r^j and r^, we find the 
connection of the Bloch wave vector Q to the total momentum of the three-particle system, P = 37Q. Using the 
periodicity of the potential Vij{r[j -t-n^jL) = T^ij(r';^), the three-particle Lippmann-Schwinger equation on a torus can 
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be written 



where the periodic three-body Green's function is given by 



or asymptotically 



1 f'^-'Jl E, 



2m J^^" E 



- q? i3 fc2 _ k?, ■ 



Next, we use the expansion of the Green's function in Eq.(B2) and the definition of the scattering amplitude in Eq.(22) 
and ( 23 1 , to obtain 



V'jAf^^(xi,X2,X3) 



2-K 



{S^JMs,r,LkMLJJM)Ys,^M',J^^,)YL'^M',^{h) 



E-m 



2m 



X i " kE,^ fs,,Lkj{kE.,) SL^Mr.^,L'^M' nL'^{kE,^rk) - Ml^ 



(Q) 



' ikE,^)jL'^{kE,,rk) 



Matching a general wavefunction of form ^^j^j cjm -0JA/(xi, X2, X3) with J2jm^JM^'jm^\^i,^2,^3), and pro- 
jecting out the partial waves, neglecting rearrangement effects from crossed channels as appropriate in an isobar 



approximation, we get a set of equations which have general form dEij F{Eij,rij,rk) = 0. Since the variables 

(fijjfk) can be chosen arbitrarily, F(Eij,rij,rk) — must be satisfied for each Eij. From this we obtain three 
determinant conditions, the first is 



det 



Ss.,Ms^^ ,S^^ A4. . cot (5s., ) - Mfl^i g, {qE^^ ) 



= 0, 



(24) 



which is Luscher's formula for scattering between z-th and j-th particles inside the isobar pair (ij). It provides the 
constraint on the phase shifts dsij as a function of invariant mass of the isobar pair {ij), Eij. The second condition. 



det 



5jM,J'M' Sl^X'^ COtSs^^L^jikE,^) 



'3 k 



= 0, 



(25) 



is a generalized Luscher's formula in moving frames for scattering between the spectator, A:-th particle, and the isobar 
{ij) with the specific spin Sij and mass Eij. Thus, Eq. ( 25 ) gives the constraint on the phase shifts ^5. j as function 
of the total energy, E, for each individual partial wave 01 isobar pair (ij), Sij. 
The final condition. 



det 



SjM,J'M' Ss^.,S'^. hk.L'^ COtds^^iqEi^) C0i5s,^Lkj{kE,j) 

- E {S[,^%;UuM^^J'M') {S,,Ms^r,L,MLjJM)Mfl,,^^^^si^,,,^ {qE.^) 



-0, 



(26) 

leads to an additional constraint on both Ss^^ and Ss^jL^J for the scattering between the fc-th particle and all the 
allowed partial waves of the isobar (ij). 

I 



In the case that only a single partial wave Sij of the isobar pair is dominant, Eq. (26) becomes redundant and 
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the three conditions reduce to two, 



cot <5s., (qb., Ms. - ,S< - Af^ . . ilE., ) 



and Eq. (25). Similar to the discussion in Section II and 



III Eq. (24 1, (251 and (261 have to be subduced accord- 



ing to irreducible representations of the appropriate little 
groups 31j. 

From Eq. (24), (25) and (26), we find that even in 
relatively simple cases, e.g. a single isobar dominating a 
single relevant partial wave, extracting the phase-shifts 
and determining the invariant mass Eij from the mea- 
sured total three-particle energy i? is a rather difficult 
task. For additional information, we can first perform 
computations of two-particle correlators with the quan- 
tum numbers of the isobar channel to obtain information 
on Eij. However, in the three-body calculation, there 
may be multiple Eij for each individual E allowed by 
kinematics (2m < Eij < E — m), thus, finding the cor- 
respondence between Eij and E may requires some as- 
sumptions and model input. 

In recent works ^4^, '5j , determining the spin of the ex- 
cited states by considering the overlap of carefully con- 
structed operators with the particular state, (njOlO), has 
been proven to be successful in certain cases. We may 
use a similar idea to identify the Eij and E relation 
in a three-particle system, by considering the overlap 
of operators with the state having the particular quan- 
tum numbers {Sij ,Lk,J) and invariant mass of the isobar 
pair Eij. For instance, working in center-of-mass frame, 
the three-particle operator may be constructed in such 
way that an isobar pair (ij) has definite relative momen- 
tum, IlPi — Pjl = ^ I rig I, rig G Z^, and definite total 
momentum, |Pi + Pj| — ^|nfc|,nfc € Z^. This opera- 
tor will strongly overlap with the state having invariant 

mass of the isobar pair (ij), Eij ~ 2^ + (^nq)^, 
and the total energy of the three-particle system E ~ 



m2-h(^nfc)2 + ^Ef^+ (^nfc)2 with smaU shifts 

caused by the interaction. Additionally, in the case of 
isobar dominance, fermion bilinear operators subduced 
from spin Sij are likely to have good overlap. 



level. Using a toy model of two-channel S'-wave scatter- 
ing we demonstrated that it is possible to determine this 
information if multiple energy levels are determined. 

Two possible strategies for extracting information from 
discrete spectra of a coupled-channel system were dis- 
cussed, one approach utilizes the near degeneracy of en- 
ergy levels in different volumes and total momenta of 
system, and another fits the discrete spectra by param- 
eterizating the scattering amplitudes with certain num- 
bers of parameters. These strategies may be useful for 
the analysis of future lattice QCD data. In particular, 
the coupled-channel analysis has to be considered for the 
strongly coupled systems, for instance, tttt, KK and rjr] 
system in S'-wave. 

For the case of three relativisitic particles undergo- 
ing scattering, we considered in finite-volume the isobar 
model approximation, where quasi-two-body scattering 
is assumed to be dominant. Under the isobar model ap- 
proximation, the three-particle partial wave scattering 
amplitudes can be factorized as the product of two indi- 
vidual scattering amplitudes, one describes the scattering 
of two particles inside isobar pair, and another describes 
the scattering between isobar pair and the spectator par- 
ticle. Three determinantal conditions, Eq. (24), (25) 



and ( 26 ) , were obtained for three-particle scattering in 
a finite volume. One condition, Eq. (24 1, relates the 



scattering phase shifts of two particles inside an isobar 
pair to the invariant mass of the isobar pair. The other 
two conditions, Eqs. (251 and (26), relate the scattering 



phase-shifts between an isobar pair and the spectator to 
the total energy of the three-particle system. A proposal 
for extracting the phase-shifts of a three-particle system 
from lattice QCD simulations is presented that makes use 
of carefully constructed two-particle operators within the 
overall three-particle operator construction. 

The finite volume representation of the isobar model 
may be suitable for systems with a sharp two-body res- 
onance. For example, in the J/ip to Stt decay [3S], the 
Dalitz distribution of the Svr shows strong bands from 
the p resonance in the sub-two-pion system, such that 
the decaying process is well described by an isobar ap- 
proximation. The rescattering effect among the three 
particles is expected to be negligible. 



VI. SUMMARY 



VII. ACKNOWLEDGMENTS 



Using the Hamlitonian formalism applied to a model 
of interacting relativistic fields, we derived a generalized 
Liischer's formula [51 [HI HH [H] for two-particle scatter- 
ing, in both the single- and coupled-channel systems, in 
moving frames. 

Our results are consistent with the ones obtained pre- 
viously in [SI [HI dH US] ■ In the coupled-channel case we 
are challenged by the fact that, even for dominance of a 
single partial-wave, the system is underconstrained for 
determination of multiple-channel phase-shifts and in- 
elasticities from a single determined finite- volume energy 



We thank D. B. Renner and M. R. Pennington for use- 
ful discussions, and our colleagues within the Hadron 
Spectrum Collaboration for their continued assistance. 
PG, RGE and JJD acknowledge support from U.S. De- 
partment of Energy contract DE-AC05-06OR23177, un- 
der which Jefferson Science Associates, LLC, manages 
and operates Jefferson Laboratory. JJD also acknowl- 
edges the support of the U.S. Department of Energy 
Early Career award contract de-sc0006765 APS ac- 
knowledges the support of the U.S. Department of En- 
ergy grant under contract DE-FG0287ER40365. 



13 



Appendix A: Relativistic Lippmann-Schwinger 
equation from Hamiltonian formalism 

Following the method presented in [2SJ [53] , we treat 
the relativistic dynamics of particle scattering in the 
Hamiltonian formalism approach. We start from the 
covariant Lagrangian in Eq.Q and choose to quantize 
the field operators in the instant form [30J ; the construc- 
tion of generators of the Poincare group can be done in 
a standard way in quantum field theory. In principle, 
one needs to solve eigenstate equations H\'if) — E\'^) 
on a instant quantization plane, where denotes the 
Poincare covariant state vector spanning the complete 
Fock space. We truncate the Fock space up to three- 
body states, assuming that this is sufficient to describe 
low-energy physics. Thus, the eigenstate equations re- 
duce to a matrix equation given in Eq.([2]). Eliminating 
the three-body sector, we end up with the relativistic 
Schrodinger-like equation for a two-body system given in 
Eq.(|3|. For simplicity, we have assumed that the two 
charged scalars scattering have equal mass, however, the 
conclusion of this work can be trivially generalized to 
non-equal masses case as well (c.f. [5^ HP] ). 

We choose the center-of-mass frame of the many-body 
system to construct multiple-particle states | JM) having 
total spin J and spin projection M. The two-particle 
state \4>'^(t>~) in the center-of-mass frame is given by 



X (27r)V!7l^(pi,P2X6;jO), 



where Pi is the momentum of the i-th particle and 
•\/s is the invariant mass of the two-particle system. 

V'jM (Pi'Pa) is the wavefunction of the two-particle sys- 
tem, describing the momentum distribution of the two 
particles. 



Similarly, the three-particle state |(/)^ 



is given by 



' ' J{27:r2Ep, (27r 



d^P2 (fp3 



^2Ep^ {2TTf2Ep 



X (27r)353(p^ +P3) 



where (/Sjjvf (Pi' P2; Ps) the wavefunction of the three- 
particle system. The wavefunctions are normalized 
so that the normalization of states is {JM\JM) = 
2y/^{2TrfS^{0). 



It is straightforward to evaluate the matrix elements 
of the eigenstate equations Eq.([2]) and to get coupled 
equations for the wavefunctions 



(27r)32v/k'2 + /i2 



y(^),(q+ik-,kO ^ y.yi,(q-ik',k- 
2v/(q + k')2 + r7^ 



(3) 



2V(q-k')^ 



2^(ik-q)2+,: 



2^(ik-Hq)2+„ 



V^-./(ik + q)2+m2-./(ik 



where we have used a short-hand notation for wavefunctions (^j^(q) and (/?j2/(q, k), with arguments of relative 
momenta defined by q = |(pi — P2),k = — P3. Eliminating the three-body wavefunction, we get a relativistic 
equation for the two-particle state \(f>~^(t>^) with an effective non-local potential generated from the neutral scalar 
exchange between two charged scalars 



with 



"(^(q.k)-'^^ ^^„^ /i,2 I ™2\ rr, ^^5— — ^ r- /^r^-. n r^. ^ n\ —. — 2' '^^^^ 



4 1 - S(q) (k2 + m2) ^(k-q)2 + ^2 _ + m2 - ^q2 ^ - ^(k-q)2+/i: 
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where the self-energy contribution is, 



S(q) 



1 



9 



1 



m? — \/o^ + m? — ^ (k' — q)2 + /i^ 



11 1 
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In coordinate space, the wave-equation becomes 



V'JM(r)= /dVGo(r-r',V^) d^'zVir' , -z) i:jM{z), 



where V'JM(r) and y(r',— z) are the Fourier transforms 
of the momentum-space two-particle wavefunction and 
effective potential, respectively. The free Green's func- 
tion is given in Eq.Q. Performing the angular integral, 
the free Green's function reads 



1 

2ir 



qdq 



-iqr 



(2^)2 ^_2^g2 + ^2' 



(A2) 




which has the following singularities in the complex q 
plane: two poles on the real axis, q = ±k, and two branch 
cuts on the imaginary axis ±[im,ioo], see Fig[7] We 
choose the contour Ci -t- C2 to include the pole at g = fc 
and circle around the cut [im, ioo] on the upper half-plane 
for first term with factor e''""'' and choose the contour 
Ci + C3 to include the pole at q = —k and circle around 
the cut [—im, —ioo] on the lower half-plane for the second 
term with factor e"*'^'", as shown in Figj?] The contour 
integral leads to 



Go(r,v^) 



1 



Pdp 

2 47rr r J (27r)2^^ 



-pr 



where k = |\/s 



4m2 is the momentum of either par- 
ticle in the rest frame of the two-particle system. The 
first term on the right hand side comes from the poles at 
q — ±fc and is proportional to the usual non-relativistic 
Green's function which oscillates over the path of propa- 
gation. The second term comes from the contribution of 
the discontinuity crossing the branch cuts at ± [im, ioo] - 
it decays exponentially over the propagation. Expanding 

^ — m^ = p^l — C'(^)^, at large separations, the 

free Green's function may be approximated by 



Go(r;ys) 



4s_ \_ 
2 47rr 



ikr 



2 e- 



(A3) 



and therefore the exponential decaying term can be 
dropped in the limit r 3> m~^ . 



FIG. 7: Integration contours and singularities of the free 
Green's function in Eq.( A2 1 on the complex q plane. 



Appendix B: Expansion of Green's function and 
regularization of expansion coefficients 

We start from the expansion of the Green's function 



1 

1? 



T - 



qePq 

k 

An 



(fcr)-^ff(3fc)j,(fcr)F,,„^.(f), 



where the summation of q runs over Pq = {q G M'^jq = 
-I- Q, for n G Z^}. The expansion coefRcients are 
given by [5] 

qePq 



fc-J q2 — fc2 



Att r 



(Bl) 



Note that the definition of nj(x) in this work differs from 
the definition in [8_ by a overall negative sign. 
Using the identities 

jXfc|r-r'|)r,„^(r^) 

^'^^ ^ ^ ^^-^'-=Jl{kr)Jl,{kr')YlraAi)Yl1m^i') 



Imi 
V mil 



.J^i±^^+^{jm,-l'mi.[lmi) {jO;l'0[lO), 
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and 

^ no(fc|r - r'l) k J2 ni (kr) {kr') Yi^a, (f ) Y^^^ (f '), 

Imi 

we also have 
1 e^q-(r-r') 

qePq ^ 

E [5i'm',Ma^ni{kr)-M^Z^_^im^k)ji{kr)_ , 

Imi 
I'mi, 

-xkj,{kr')Yir^,{v)Y*^^,{i'), (B2) 



with 

X /^'^±i|^(jm,;Z'™rKm;)(iO;Z'0|ZO). (B3) 

If Q is identified with for two-particle scattering 
and for three-particle scattering, these expressions 
are the same as those presented in [9 , and the regulari- 
sation procedure outlined therein can be followed. 
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